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VI. The Scattering of Plane Electric Waves by Spheres, 
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Introductory Note. 

The problem which gives its title to the present paper has been handled by various 
writers, notably by Lord Rayleigh, Sir J. J. Thomson, and Prof. Love. In most 
cases the solutions have been expressed in a Cartesian form ; but it appears to me 
that a marked simplification is introduced by using spherical polar co-ordinates. The 
preliminary analysis becomes shorter, and the conclusions are easier to interpret ; in 
fact, the analysis is nearly as simple as in the analogous problem of electrostatics, 
when an electric field is disturbed by the presence of a dielectric sphere. 

To obtain the requisite solutions a new general solution of the electromagnetic 
equations in Cartesian form is given in § 1, and is then transformed to the spherical 
polar form ; §§2,3 contain a summary of certain analytical results required in the 
sequel. 

§ 4 contains the general solution of the problem of finding the scattered waves 
when a plane simple harmonic wave strikes a sphere ; and in § 5 the solution is 
applied to the case of a small sphere. These formulae (all of § 4 and part of § 5) were 
originally worked out in 1899, but publication was postponed in the hope of 
completing the problem of the large sphere. 

In § 6 the problem of a large sphere is considered by applying to the formulae of § 4 
a method of approximation devised by Prof. H. M. Macdonald # for dealing with 
waves incident from a Hertzian oscillator on a conducting sphere. The formulae of 
§6 were worked out early in 1910 and were given in my University lectures at 
Cambridge in that year . t 

At the same time I succeeded in obtaining a different treatment (given in § 7 
below) which confirmed the other results, and gave an easier process for dealing with 

* 'Phil. Trans. Roy. Soe.,' A, vol. 210, 1910, p. 113. Prof. Macdonald tells me that he had worked 
out (at about the same time) results in reference to the problem of §6; but these have not been 
published. 

t An alternative solution was obtained by Prof. J. W. Nicholson at about the same time; his solution 
starts from Sir J. J. Thomson's formulae. Prof. Nicholson's results originally differed from those of 
<§ 6 ; but on revision agreement was obtained (* Proc. Lond. Math. Soc.,' vol. 9, 1910, p. 67 ; vol. 11, 1912, 
p. 277). 
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points behind the sphere. The method of this section is similar in some respects to 
one used by Prof. Macdonald in a later paper. # 

The formulae of §§6, 7 have been delayed in publication for two reasons : in the 
first place I wished to obtain some confirmation from direct numerical calculation. 
This has now been carried out by Messrs. Proudman, Doodson and Kennedy (of 
Liverpool University).! It appears that the agreement with the formulae of § 6 is 
quite close (for kco = 9, 10) from 6 = 0° to 90°, and for the Z-component up to about 
120°. The formulae of § 7 also give good results in a cone of about 10° behind the 
sphere (that is, from 6 = 170° to 180°). It is clear, however, that an approximation 
suitable from 6 — 90° to 170° (for Y) and from 6 = 120° to 170° (for Z) has still to be 
obtained. But nevertheless the present approximations proved a valuable auxiliary^ 
in checking and testing the numerical work. 

[The paper in its original form was presented to the Society on April 13, 1916 ; 
owing to the difficulties in regard to labour and paper during the war, I was asked to 
condense the introductory matter of §§ 1—3. This proved to be impossible until now, 
on account of pressure of war-work of various kinds. In the present version § 1 has 
been re- written so as to reduce its bulk ; in §§ 2, 3 certain formulae have been omitted 
which were not used in the applications of §§ 4-6. 

In re-arranging the paper it proved convenient also to number the formulae 
differently. The decimal system has now been adopted ; here the figure before the 
decimal point indicates the section of the paper in which the formula occurs. The 
figures following the decimal point are to be regarded as following the same order as 
ordinary decimal fractions. Thus (5*21) and (5*22) fall between (5*2) and (5*3), and 
all these formulae occur in § 5. — Added March 18, 1919.] 

§ 1 . A General Solution of the Fundamental Electromagnetic 

Equations.^ 



The fundamental equations of electromagnetic waves may be written 
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* 'Phil. Trans. Roy. Soc./ A, vol. 212, 1912, p. 299. The two methods are not identical; but they 
appear to yield equivalent results in all the cases to which they have been applied. 

t 'Phil. Trans. Eoy. Soc./ A, vol. 217, 1917, p. 279. The calculation was originally undertaken by 
Dr. Proudman in consequence of a suggestion made in my lectures of 1912; the work, however, proved 
to be longer than had been anticipated and was completed by Messrs. Doodson and Kennedy. 

J See the paper last quoted, p. 292 et seq. 

§ Revised March 18, 1919 ; see note at the end of the introductory remarks above. 
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Here (X, Y, Z) denotes the electric force, (a, /3, y) the magnetic force, K is the 
dielectric constant, /x is the magnetic permeability, and the axes of reference are a 
Cartesian right-handed system. The units adopted are those of the electromagnetic 
system, and c is the fundamental constant generally identified with the velocity of 
radiation in free space ; the equations (E) are those derived from Ampere's law, and 
the equations (M) are similarly derived from Faraday's law, the two together 
constituting the circuital relations of the electromagnetic field. 

It has proved possible to obtain a solution of a very general type, by assuming that 

(i*i) a — x ®' a — y ^ a — 

then equations (M) yield 

{l2) f "dt~ y dz "dy "dt~ Z dx X dz M 3i X dy V dx 

Substitute from equations (1*2) in the first equation (E) and we obtain 

pi2 £\2 £\2 

where A 2 denotes Laplace's operator ^ + ^ + ^-§ • 

On comparing equations (1*1) and (l'3) they will be seen to be consistent provided 
that 

^*> c 2 df ~^ + X dx +y dy +Z dz 

and that 

(1-5) ^|S = A 2 Q. 

Thus Q must satisfy the fundamental wave- equation, which is satisfied by any 
component of the electric or magnetic forces (X, Y, Z) or (a, /3, y). 

For our purpose it is more convenient to express the above solutions in terms of 
spherical polar co-ordinates r, 6, cj> ; these are supposed to form a right-handed system, 
when taken in this order, so as to avoid changes of sign in introducing the new 
co-ordinates. We write here # (R l9 K 2 , E 3 ) for the components of electric force in the 
directions of r, 6, cj> respectively ; and (H 1? H 2 , H 3 ) for the components of magnetic 
force. 

Equations (l'l) then become 

* This is done to avoid confusion with the Cartesian components used in equations (E) and (M) ; but 
in the subsequent sections we shall use (X, Y, Z) and (a, )3, y) for the spherical polar components 
here denoted by (Ri, R2, R3) and (Hi, H 2 , H 3 ) respectively. 

2 c 2 
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As regards the transformation of (l'l) to (I'll), it is sufficient to note that the 
gradient of P has the spherical polar components 

3P iap 1 3P 



•? — 



dr r d6 r sin 6 d<j> 

and that (r, 0, 0) corresponds to the Cartesian vector (x, y, %). 

To obtain the formulse corresponding to (1*2) we observe that the vector on the 
right is equal to the vector-product of the two vectors, 

(*, y, z) and (g, |Q, |2); 

and that these two are represented by 

(r, 0, 0) and (f, 1 «jl -i- ® 

\ or r ad r sm 6 0$ 

Thus the vector-product has the spherical polar components 

0, _^L_9Q, 3Q 1 

sin 6 d<j> dd. 

Consequently equations (1*2) now become 



(1'21) - M ^il=0, 



3 Hi 3H 2 1 8Q 3H 3 3Q 

a* ' M "3« sines/ M a« 3# 



while (1*4) and (l"5) give 

(vri\ M K3 2 Q_ 1 3 / 2 3Q\ 1 3 / ■ fl 3QV 1 8 2 Q 

1 ' V 3£ 2 ~ r 2 drY dr) + r 2 sin ddY™ 30/ r 2 sin 2 3^> 2 

-lil/n^ 1 ' 3 / • fl 3Q\ 1 3 2 Q 

- r 3r * ^ + r 2 sin 30 ^ sm y 30 J + ^ arfe fy» • 

A consideration of these formula? suggests that further simplifications can be 
obtained by writing 

(re) p-£, Hj-f 3 ^. 

which together satisfy equation (l'4l); and then equation (l'5l) leads to the 
equation for U : — 

(V7\ mK 3 2 U _ 3 2 U- 1 3 / • 3U\ 1 3 2 U 

{ } c 2 dt 2 ~ dr 2 + r 2 sin ddV™ dd J r 2 sin 2 6 d<p 2 
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Substituting from (l'6) in (I'll) and ( 1*21) we obtain the final expression for the 
field in terms of U :— 



r 



(1-8) 
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_ a 2 u m k a 2 u 



dr 2 

i a 2 u 
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r d6 dr 

i a 2 u 

p d<j> dr 



cHi = 

tt - ^lA/Kau' 

C±±2 ~ + P d<p\c dt 



p = r sin 0, 



cH, = — 



1_3 /K3U N 



In like manner we obtain another set of solutions by making an assumption similar to 
(ri) for the components of magnetic force (a, /3 9 y). This gives the field : — 



f 



(1-9) 



R 1 = 



*S 



Ro = 



R» = 



_ 1 d_U ST 

p d<j>\c dt 

1_3_/m3V n 
rdd\c dt 



cHj = 



cH, = 



cHo = - 



a 2 v m k a s v 






ar 2 c 2 a* 2 






i a 2 v 

r dOdr 


P = 


r sin 0, 


i a s v 






p d<pdr 







where V is a second solution of equation (l"7). 

It can be proved* ..that "(1*8) gives the most general field in which the radial 
magnetic force (E^) is zero, while (l"9) gives the most general field in which the 
radial electric force (R x ) is zero. It can also be shown that the field is uniquely 
determined by the value of R x and H^ ; and accordingly the most general solution 
can be obtained by the superposition of (l"8) and (1*9). 



§2. Further Specialization of the Solution of § 1. 

If we superpose the fields (V8)< (l'9), and now utilize (X, Y, Z), (a, ft y) to denote 
the spherical polar components of the field, we have the general solution f:— 



r 



(21) 



x = a 2 U_^K3 2 TJ 



dr 



dt' 



■ Y = i a 2 u i a (fj.dv 



r dddr p d<p \c dt 



p = r sin 0, 



Z = - 



i a 2 u . i _^ a* av\ 



+ 



pd<j>dr r d0 \c dt 

* See a paper in the 'Philosophical Magazine/ July, 1919 (6th ser., vol. 38), p. 143. 
t Originally worked out in 1899, and first published as a question in Part II. of the ' Mathematical 
Tripos/ 1910. 
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3 2 V fxK 3 2 V 

Co. = 



(2-2) 



9r 2 c 2 dt 2 



. 1 3 2 V , 1 3 /K3U\ 



\. 



c 7 



r d6 dr p d<f>\c dt / ? 

1 d 2 Y 13 /K 3U\ 7 
p 3c/> dr r d6\c dt J 



where IT, V are any two solutions of the equation 



p = r sin 0, 



(2*3) 



M K 3 2 U 3 2 U 



+ 



( Sill -^-r I + 






3 2 U 



c 2 dt 2 dr 2 r 2 sin 6 dd\ ^ dO ] r 2 sin 2 6 d(p 



A solution of (2*3) which is sufficiently general for the applications in view may be 
found by assuming that U and V can be expressed as sums of terms of the type 

F(r, t)xY(0, 0). 

It is easy to see that then (2'3) leads to the equation 



(2'4) 



F\3r 2 c 2 dt 2 , 






( sm 9 -z-r + 



1 3 2 Y~ 



Y I sin 6 dd \ 30 / sin 2 d dtf 



and since the two sides of equation (2'4) are functions of r 9 t and of Q, <j> 
respectively, it is clear that each side must be a mere constant. If we write the 
constant in the form n(n+l), it is evident that Y must be a surface-harmonic of 
order n. 

Accordingly in problems (such as those with which we shall be concerned in the 
sequel) where the whole of angular space is considered, the value of n must be a 
positive integer ; for (except when n is an integer) there are no surface-harmonics 
which are everywhere continuous and single- valued. 

Thus we may reduce our solution to the form 



(2*5) 



IT or V = 2F n (r ? t)Y n (6 9 <p), n = 0, 1, 2, 3, ..., 



where F n is a solution of the equation 



(2-6) 



and 



3 2 F. 



n 



dr" 



1 3 2 F n n(n+l) ^ __ 

2 ^j.2 2 n ~~~~ ' 

cr at r* 



c 2 = c 2 /( M K). 
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The general solution of equation (2 '6) is well known, and it is given by* 
(27) F n (r, «) = r- 1 (- \ A)' j /^-^+y^ + r) ' 

where the functions / and gr are arbitrary. 

In the special case of divergent waves, the function g can be omitted in (2*7) ; and, 
if the region considered includes the origin, then gicj + r) = —/(c^ + r), so as to 
make F n (r, t) continuous at r = 0. 

It will be convenient to notice that in consequence of equation (2*6) the radial 
components of force can be written in the simpler forms 

(2-8) X or « = 2^F.(r I i)Y,(() ) 4 

It will be noticed that we can at once determine the form (2*5) for U or V 
when the radial forces have been expressed in the form (2*8) ; this agrees with the 
general conclusion stated at the end of § 1, that (in spherical polar co-ordinates) the 
remaining components of force are completely determined when the two radial 
components are known. 

§ 3. Special Case of Simple Harmonic Waves and the Appropriate 

Functions. 

We assume in future that the waves are simple harmonic, of wave-length 2-7t/k in 
free space ; we can then suppose the time to occur only in the form of a time-factor 
e mc \ with the usual convention that finally only the real (or the imaginary) parts of 
the formulae will be used. 

The functions f, g occurring in equation (27) above are then exponentials of the 
types 

e «i<**-r> an( J giiito+r)^ 

where at 3 is given by 

/CjCi — /cC, Or /q = k y^/xK). 

Thus (if we now suppress the time-factor e tKct ) the functions given by (27) are 
of the types 

(3-i) r „ +1 / i aye-*" rM+1 / i a_ygn r , 

\ v drj r \ v dr/ r 

We shall be concerned with two special types only : (i.) divergent waves ; 
(ii.) waves which are continuous at r = 0. The former of these corresponds to the 

* See, for instance, Lamb's ' Hydrodynamics/ 1906, art. 295; an alternative method of solution 
given in § 3 of my paper in the ' Philosophical Magazine/ quoted on p. 179 above ; compare also A. E. H. 
Love (' Phil, Trans. Roy. Soc./ A, vol. 197, 1901, pp. 9, 10). 
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first expression in (3'l) ; while the latter is found by combining the two expressions 
so as to yield 

(3-2) rW-1 J-Y ^isd . 

\ r or J r 

Using the notation explained in (3*4) and (3*5) below, the standard functions are 

E n (/^r) for divergent waves, 

ano 

S n (^r) for waves within a spherical boundary. 

Consequently, for waves inside a spherical boundary, (2' 5) and (2'8) can now be 
replaced by the forms 



(3*3) < 



IT or V = 2S n ( Kl r) Y n (0, <p\ 

X or ca = 2 — ^—g — - S n (/qr) Y tt (0, <f>) 



for divergent waves the function S n (/c x r) must be replaced by E n (/qr). 

Definitions and Properties of the Two Standard Functions S w (z), E n (z). 
We write for brevity 

(3'4) S n (z) = z n+lf l d \ n /* m * 



'4, — * 



£ <W \ £ 



3* +l L 3* S* 






1.3. 5„.(2w + l) I 2(2^ + 3) 2. 4 (2^ + 3) (2^ + 5) "j* 
In terms of the known Bessel function we can write 

<»■«* s.W = V(f)j-W. 

and accordingly the fuaction S ?l (z) is the same as that denoted by u in one of 
Macdonald's papers. # 

In the notation adopted by LAMB,t and those writers who have used Lamb's 
solutions as the fundamental forms, we have the identity 

(3*42) S n (z) = z^ n {z). 

* 'Phil. Trans. Koy. Soc./ A, vol. 210, 1910, p. 113. See in particular p. 115, 
f < Hydrodynamics,' 1906, Art. 287. 
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Similarly, we write 

(3-a) E-M^'B !)"(?' 

= C n (s)— £.S M (s), 
where 

(3-6) GJz) = z^t- 1 - #V (™** 

\ z dzj \ z 

_ 1 . 3...(2w-l) L z 2 z* 1 

z n 1 2(l.-2w) 2.4(l-2w)(3-2n) "T 

In terms of the K„ function (the modified Bessel function used by Maodonald) 
we have the relation 

(3-51) E„ (z) = V(?) e * M)mK "+i (")• 

Thus 

E w (z) = v— m 

in terms of the notation used by Maodonald in the paper last quoted, and in 
Lamb's notation 

(37) E n (z) = ^ +1 {* n (s)-!^ W} = ^ +1 /n (*)• 

In consequence of the equation (2*6) we see that both S n (z) and E n (f) are solutions 
of the differential equation 

(3-8) ^ + {i_?L(?L+Hjs n =0, 

The functions S n (s), C n (s) and |E n (s)| have been tabulated from z = 1 to 10, and 
for values of n ranging from to 22, by Mr. Doodson, # and these tables have formed 
the basis of the numerical calculations mentioned on p. 176 above. f 

It will be convenient to collect here the simple relations amongst the functions 
S n-1 , S n , S n+1 , which correspond to the known results for Bessel functions, or to those 
given by Lamb for the equivalent function \fs n (z). 

Difference Relations for the Functions S n , E n . 
From (3 ; 4) we see that 



"^n 



dz 



(3-8i) s„ +l( .) - _*«.jjL{§k&} = •lis.w- 

and by using (3 • 5) we see that the same relation holds for E n (z). 
Again, it will be found that 

[§L + 5Wp « (2n + v) z n+ P~\ 
\dz zj 

* 'British Association Beport,' 1914. 

t Proudman, Doodson and Kennedy, 'Phil. Trans. Eoy. Soc./ A, vol. 217, 1917, p. 279. 
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and using this in equation (3 * 4) we deduce that 

(3 • 82) S„_! (*) = (| + j) S M (z) = I S B (z) + ^ ■ 

Combining (3-81) and (3 -82) we have also 
(3 • 83) S ?l _! (z) + S n+1 (z) = 2n+1 S n («). 

Z 

The relations (3-82), (3*83) hold equally for E n (z) and C n (z), as may be seen from (3 -5) and (3*6). 
As E n (#), 8 n (z) are independent solutions of the equation (3 * 8), it is evident that 

E H (z) — ~ - S n (*) ^2 = const. 

dz dz 

Now when z is small, it is easy to verify from (3 -4) and (3*6) that 

E " W f - Si + ° ^ s » W f" - " 2ST1 + ° <*>' 

and accordingly we have 

(3 • 84) E„ («) ^ - S„ (2) ^j* = 1. ' 

In the discussions of § 6, when n, z are both large, it will be convenient to adopt 
the following notation : — - 

(3*85) | E n (z) | = E, E, (z) = He- 1 *, so that S B (0) = E sin ^, C n (2) - E cos ^. 
Substituting from (3*85) in (3*84) we deduce that 

(3*86) K*^£ = 1. 

dz 

Before leaving these preliminary formulae it will be convenient to quote the formula 
for e lKZ in terms of our standard functions ; namely 



(3*9) 6"*= 2 (2w + l)i»^^P n ( M ), 

where z = r cos = 77* and P n ( M ) is Legendre's polynomial of order n. 

This result follows at once from the formula given in Lamb's ' Hydrodynamics/ 
Art. 291, on using the relation (3*4) between ^(kt) and S B (*r), already quoted. 
It is of course evident that an expansion of the type (3*9) might be anticipated, 
since each side satisfies the wave-equation, is symmetrical about the axis of z, and is 
continuous at r = ; the determination of the numerical coefficients may be then 
carried out quickly by comparing the terms in (*?>)* on the two sides of the equation. 

§4. Plane Electromagnetic Waves Incident on a Spherical Obstacle. 

Suppose that the incident wave-train is travelling along the negative direction of 
the axis of z (that is, from = towards 6 = tt) ; and that it is polarized in the 
plane of yz (that is, in the plane <p = |tt). Suppose further that the electric force in 
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the wave-train has unit amplitude ; then, in terms of the Cartesian specification, 
the incident wave is defined by # 

X = - e lK(ct+z \ cp = + e LK{ct+z \ 

the remaining components of force being zero. 

We must first express this wave in the standard forms of (2*1) and (2*2); we 
therefore introduce polar co-ordinates, and then proceed to find the radial components 
of force, which will suffice to determine the functions IT, V. 

These radial components are given by 



(4-1) 



X = —sin 6 cos (j> e 



IKT cos 9 



ca, = + sin 6 sin <f> e 



iKf cos Q 



where the time-factor e lKCt is now omitted. 
Now from (3*9) we have the formula 



e 



<-*?• cos Q 



= £ (2n+l)i n ^^P n (coaO). 

71 = KV 



So, differentiating with regard to 0, we find that 



(4-2) 



singe 1 '"' 0039 



OD 



2 (2w+ 1) L n ~ l Mtt P'» ( cos e ) sin e - 



Accordingly, on substituting (4*2) in (4*1), we find that in the incident wave 



r 



(4-3) ■< 



and 



u = _sin0coB_^ - (^±li t »-iS,(^)F fl (co8 0) > 

r n = in(n+l) 



V 



K 



n 



in{n-\-l) 



by comparing the two formulae (2*5) and (2*8). :: 

The corresponding waves in the interior of the sphere will be given by the two 
functions 



r 



(4-4) 



u Bine cob* £ ^il r!BAMF>(cos0) , 



K 



n 



= i n 



< 



and 



(n+1) 



Vi = 



V. 



where 



+ sin0cos^ - 2n+l ^^^ )p , (oos0)! 

AC 2 n = lW(W+l) 



and K, //, are the fundamental constants of the spherical obstacle. 

* It is assumed that in the incident wave we may take //. = 1, K = 1, £i = £5 *i = *< 

2 D 2 
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Similarly the scattered waves will be given by the two functions 



r 



(4*5) \ and 



TT sin 9 cos <f> * 2n + l n ~i a t? / \ r>/ / m 

(j o — _ _ r 2, — -? ~v * A n hi ft \kt) r n (cos 0) 

/c n = i n{n-t-l) 













sin sin <j> 



oo 



k 



n 



i - 1 W 



(n + l) 

2n+ 1 
(^4-1) 



to- 



" 1 C n E n (kt) P' n (cos 0). 



The boundary conditions are given by the continuity of the tangential components 
of electric and magnetic force at the sphere r = a. 

It is evident from the form of equations (2"l), (2*2) that these conditions will be 
satisfied if we take 



r 



(4-6) { and 



V. 



8U, = 9U 8U n 

dr dr dr 



8V, = 8V 8V» 
dr dr dr 



KU, = U+U c 



fi 



Vi = v+v . 



Thus we find that A n and C n (the coefficients in the scattered waves) are given by 



r 



(47) \ and 



{S n (kg) -f A n E n ( K a)} -% ■ n ) K ia { = S' n (/ca) + A B E' n (*a) 

/civ D n \^i^/ 

/C// O n (/Cj(Xj 



The special case of a perfectly conducting sphere is given by making the tangential 
electric force zero at the sphere r = a ; and this condition is satisfied if 



(4*61) 






- v+v, 



Thus we find the simpler formulae 



(471) 



S' B (tea) + A W E' W ( K a) = S n (*a) + C W E„ (*a) = 0. 



which may be regarded as limiting forms of (47), when | K | -> go and /a -> 0. 

The formulae (4*5) and (471 ) lead at once to those quoted by Dr. J. Proudman* 
in calculating the pressure of radiation due to a plane wave incident on a small 
conducting sphere. 

In all the applications with which we shall be concerned at present the point at 
which the disturbance is to be calculated will be at a distance large compared with 



* ( 



Monthly Notices of the Royal Astronomical Society/ vol. 73, 1913, p. 535, 
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the wave-length- Then we can simplify the general formulae by observing that (3*5) 
may be replaced by the approximation 



,— tur 



E n (kt) = i n e- 

if l/xr, (1/kt)'\ &c., are neglected. Further, in the final formulae for the forces, 
U and V occur only in the two combinations 

-» jr ___ 1 8U = ikU -^ _ l_^Vo __ 

n* 7)^* ty* r or r 

where terms of the relative order I/kt have been rejected. 

On substituting from (4*5) we find, to the same degree of accuracy, 



r 



(4-8) 



< 



V 



M= +i^ = --U o = sinflcos0— 2 (-i)» 2 /* + j ; A n F n ( cos 0) 
r dr r kt n = i n\n+l) 

and 

N=-i^ = +^V o = Bin0Bin0^^ I (-1)- 2n+l CF (oos ^ 



Then, substituting in the general formulae (2*1) and (2*2), we find that (to our 
present order of approximation) the radial components of force are zero, and that the 
transverse components are given by 



(4*9) i 



V 



v- 3M 1 3N 

30 sin 6 defy Y 

rj 1 3M- 3N Q 

sin 6 d(p 66 



Accordingly the electric and magnetic forces in the scattered waves are at right 
angles to each other and to the radius, and their magnitudes are related in the same 
manner as in a plane wave. 

This conclusion might very well have been anticipated ; and for the case of small 
obstacles of any shape (with constants K, ju differing but little from unity) the 
conclusion is contained in a paper by Lord Rayleigh. # But I cannot find that it has 
been noticed for the case of spheres of any size, and of any electrical and magnetic 

constants. 

This may serve to indicate one advantage of the formulae in spherical polars over 

those in Cartesian co-ordinates. 

The formulae (4*8), (4*9), with the values of A H , C w given by (471), were those 
used by Messrs. Proudman, Doodson, and Kennedy in their numerical calculations 
quoted in the introduction to this paper. 

* ' Scientific Papers/ vol. 1, pp. 522-536. For a small perfectly conducting sphere the same conclusion 
is given by Sir J. J, Thomson, < Recent Researches/ p. 448, 
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§ 5. Spheres Small Compared with the Waye-Length. 

The fundamental assumption is that | tea | is small enough to justify us in rejecting 
all but one or two terms in the power-series for S n (ko) and E n (ko). 

It has been usual to assume further that | K X a \ is correspondingly small ; but 
Dr. Proudman has remarked that in the case of a dielectric sphere with a large 
value of K the second assumption need not follow from the first. It seems worth 
while therefore to simplify the formulae (4*7) by expanding in powers of |kc&|, while 
retaining the general forms for S n (^a) ; it will be seen moreover that the resulting 
formulas form a link between the results for dielectric spheres and those for 
conductors. 

Remembering that S n (/ca) is of order (/ca) n+1 and that E n (/ca) is of order (/ca)~ n , it is 
easy to see from (4*7) that in general A n and G n are both of order (tca) 2n+l . Thus in 
the first approximation it will be sufficient to deal only with the coefficients Ai and 
C x ; and for these we need the formulae for S x (^a) and S'i (/qa). Now from (3*4) 
we have 



L d 



ana so 



Hence 



where now 



a / \ _ 2/ _ I d_\ /sin z\ _ sins __ __ sins ,- x 

KJi \as) — <0 i ^r^ I I / — ■ — OUd as — \ X 6 bUl As /, 

' x ' z dzj \ z I z z 



el/ / \ • ( .. 1 \ , cos z sin z ( 2 ' , . \ 
b i (2;) = sin z( 1 -j 4 = — j— \zr—l+z cot 0). 

Z I Z As 



zS ; 3 (z) ___ z 2 —l + z cot Z __ -p / \ ___ . 

Si (z) 1—zcotz 



F(z) = z 2 /(l-zcotz). 



In the second place, for the functions of /ca, from (3*4) and (3*6) we find the first 
approximations 

S 1 («) = K E 1 («) = l/s. 

Substituting, it will be seen that for n — 1 the first equation in (4*7) gives 



l(/ca) 3 j K (ko,) s 

where 

^M = F ( tia ). L 

bi \K\a) 
After a little reduction the last equation gives 

The second equation in (4*7) gives a similar formula for C 1? with /* taking the 
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place of K. However, in most cases, it is sufficient to write jm = 1, and then the 
formula simplifies further and becomes 

(5-2) 0, - M. {^-j -i 

The two formula (5*1) and (5'2) are due to Dr. Proudman, who has pointed 
out that they connect the results found by Lord Rayleigh for the case of dielectric 
spheres, and by Sir J. J. Thomson for conducting spheres. 

To deal with the ease of dielectric spheres we do not regard K as large, so that K^a may be regarded as 
small (of the same order as ko) ; and then the approximations 

1 - (*ci&) cot (kjo) = |('<i#) 2 , or F(/qa) = 3, 

may be used. This gives, in place. of (5*1), (5*2) the simpler forms due to Lord Rayleigh* 

(5-21) A a = §(*<*)» |^I, d = 0. 

On the other hand, Sir J. J. Thomson's case corresponds to the assumption that K is of the form 
Ki - iK 2 where K 9l is very large ; then | K X a | may be regarded as large, and K Y a as complex, with a 
negative imaginary part. Thus approximately! cot (/q<z) = i, and so \¥ (k^)] = | /c 1 «, | , which (although 
large) is small compared with | K | = | K^a \ 2 /(ko) 2 . Hence we find from (5 * 1 ) and (5 • 2) the approximate 
results, 

(5-22) A a = §(**)», d = -|M 3 

as given by Sir J. J. Thomson. J Of course this pair of formulae follow at once from (4*71), on inserting 
the approximations for Si (ko) and Ei (ko) given on p. 188 above. 

Dr. Proudmajst makes the further remark that, under the conditions assumed in 
(5*1) and (5*2), variations in the wave-length may produce very considerable changes 
in the magnitudes of A x and Ci, on account of the presence in F (^a) of cot (^a)? 
which may vary very fast. It is of course supposed that the sphere is dielectric, 
otherwise cot(/c a a) could be replaced by <, as already stated, t 

It is worth while to note the simple formulae for the scattered wave, derived from 
(4*9) ; these give, to the present order of approximation 



f p-lKf 



(5*3) 



Y = +Cy — (fC^— I A! cos 6) cos <j>, 



e- iKr 



Z = —c8 = (f A 1 — -|Ci cos 6) sin <b ; 

kT 



* * Scientific Papers, 7 vol. 4, p. 321 (106) ; see also vol. 1, p. 526. 

t Provided that the imaginary part of /q<2 exceeds ir in numerical value, the error in this approxima- 
tion is less than one half per cent. 
\ < Recent Researches,' p. 448. 
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and so (to this order) the scattered wave is zero in the direction given by 
(5*31) = 0, coserrCi/Ax, 

provided that A x is numerically greater than P Thus in Lord Rayleigh's case, the 
direction is given by 9 = \tv ; # and in Sir J. J. Thomson's by — § tt. 

It may be noted here that, if the sphere has a sufficiently large dielectric constant K, it may happen 
that Ai is numerically less than Cj ; and then the direction given by (5*31) is no longer real. 

Taking K to be real (the case of a conductor having been already considered on p. 189), it is easy to see 
that A l < C T gives F (*-,&) < 1 (on the assumption K > 1). Now the function F (z) steadily decreases 
from 3 to as z varies from to w ; and a rough calculation shows that for z = |^r, F (z) is slightly less 
than unity. Also in order to justify the approximations used for S 2 (ko) and E x (ko), we must suppose 
that (ko) 2 ^ ^q. 

Hence the possibility contemplated may occur if, say, 

(|tt) 2 < ( Kl af < 77-2, and (mf S T V> 
giving 

K>75. 

The direction in which the scattered wave vanishes will be given by 

(5 -32) 4> - |tt, cos = A^ - 2F (/c^)/{3 - F (^a)}, 

the final formula being simplified by remembering that K is large. 

I am not aware that there is any experimental evidence showing traces of this phenomenon ; in fact all 
the evidence shows that <£ = 0, = ^w is not far from the truth. Thus the circumstances in actual 
experiments cannot have been such as to introduce the reversal of magnitude between A Y and G v 

(ii.) Second Approximations. 

We proceed next to find second approximations, assuming that | K | is not large ; 
it will be necessary to retain the second terms in the series (3*4) and (3*6) for S x and 
E l5 but the first terms will suffice for S 2 and E 2 .f It is easy to see that then terms 
of order (/ca) 5 occur in the coefficients A l5 C x and A. J? but that no other coefficients can 
contain terms of order lower than (/ca) 7 . 

Using now the series (3*4) for S x (z) we have 

s x (z) = ^(i-M s\ (z) = §* (i -&), 

retaining the second terms only in each series. Thus we find, to the same order, 

b i (z) i 2 

h x {z) 
Also (3*6) gives similarly 

z z 

* A closer approximation is worked out on the next page; see formulae (5*6) below, 
t It would be possible, of course, to obtain second approximations to (5*1) and (5*2), but a glance at 
the formulae shows that the work is so laborious as to be almost impracticable. 
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On substituting these results (47) becomes 



K 5 



and 



2- 



Kp : 



fM+f™ 



ii 1 -®^ 1 ** 



3 V 5 



p\ 2 



2^ 
3 



^V^i-f 



where, for brevity, we have written 



ic« = p, (^i^) 2 = Kp s 



On reducing these equations, the results are 



r 



K-l 



(5'4) < 



xV ~~~ Zi \ 9 

p 2 



and 



Ai = l ^frtM 1+l \K + 2 



CV=A(K-l)i> 6 . 



To determine A 2 , the series for S 2 (z) and E 2 (2) will be required to the first terms 
only ; these are 

S 2^) = tz, E 2 (z) = -, 



15 



z 



and on substituting in (4* 7), we find 



(5-5) 



2 15V2K + 3, 



The field of the scattered waves is then given by 



(5-6) 



e 



— LKV 



Y = +cy = (fCi-f A x cos 0+f A 2 cos 20) cos 5 



^ 



Z = -c/3 



e' 



•IKf* 



V 



KT 



(f A x — f C x cos 0— f A 2 cos 0) sin <£. 



The field (5*6) is accordingly zero (to the same degree of approximation) in the 
direction given by 

# = 0, cos = (CWAO/A, = A (K ~^^ +2) M 2 - 

This conclusion is apparently new; but it confirms an approximate result due to 
Lord Eayleigh, # according to which the scattered wave is zero in the direction 

given by 

= 0, cos 6 = 2 X 5 (K-l) (/ca) 2 , 

w&en (K-l) is treated as small. But on the other hand, our result contradicts 

* 'Scientific Papers/ vol. 1, p. 531, formula (61). 
VOL. CCXX. A. 2 ' E 
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a statement made by Prof. Love # that there is no direction in which the scattered 
wave is completely cut out ; however, on a closer examination of Prof. Love's formulae, 
they appear to confirm the present conclusion. 

The formulae in question are (42) and (43) of the paper just quoted, but apparently there is a slip in (43). 
In the last line of (43) the factor given as should really be (z 2 ~ y 2 )/r 2 ; the source of the inaccuracy 
being apparently in the passage from the formula (39) to (41). On introducing this additional term in 
the magnetic force, it appears that the electric and magnetic forces are zero in the direction given byt 

* -U ' r ~ 15(2K + 3) {a) 
in Prof. Love's notation ; of course " zero " means that the forces are really of order (*&) 7 at most. 

It is not difficult to prove that the formulae (5*6) agree with those found by Lord 
RayleighJ and Prof. Love ;§ the method to be adopted is similar to that used in § 4 
of my paper in the ( Philosophical Magazine ' (quoted on p. 179 above). But it should 
be observed that in the specification of the incident wave adopted by Lord Rayleigh 
and Prof. Love, the electric force is parallel to the axis of y ; but here the electric 
force is parallel to the negative direction of x. Thus if <$> denotes the azimuthal angle 
corresponding to the former specification, it is evident that $ = \ir corresponds to 
(j> = it ; and accordingly we shall have in general the relation 

because both angles are measured in the right-handed sense about the axis of z. 

Lord Rayleigh's paper contains tables and graphs from which it is easy to 
determine the variation of the field with 6 ; and in order to connect his tables with 
our formulae, let us write (5 '6) in the form 



r eT lKr e~ lKT 



(57) \ 



\. 



Y = +cy = R cos 0, Z = — c/3 — S sin <j> 9 

kT kT 

where 

R = f d-fAj cos + £ A 2 cos 20, 

S = f A ± - f d cos 6 - f A 2 cos e. 



Consider, first, points in the plane given by x = 0, in Lord Rayleigh's notation ; 
this gives <p f = Jtt, or d> = tt. Hence, omitting the factor e~ lKr /(Kr), the electric force 
is equal to R (in the direction of 6 decreasing) ; and accordingly R is represented by 
the graph of the Cartesian component (yZ—zY)fr, given by Lord Rayleigh. 

Secondly, consider the plane y = ; that is <j> = 0, or <j> = |-tt. Here the electric 

* 'Proc. Lond. Math. Soc./ vol. 30, 1899, p. 318. 

t A numerical slip in the first line of each of the formulae (42) and (43) has to be corrected; the 
correction was given in the Errata, vol. 31, * Proc. Lond. Math. Soc.' 
J l Scientific Papers/ vol. 5, p. 559, (u) and (v). 
§ Formulge (42) and (43) of the paper just quoted (allowing for the corrections just mentioned). 
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force is equal to S, perpendicular to the plane ; thus S is represented by the graph of 
the Cartesian component Y. 

In general, the resultant electric force is represented by 



e" lKT 



Kr 



\/(R? cos 2 <p + S 2 sin 2 <£). 



§ 5. (iii.) Second Approximations for Conducting Spheres. 

The foregoing algebra needs no alteration beyond replacing y^K by the appropriate 
complex refractive index associated with the particular metal and wave-length 
considered. This of course assumes that | K | is not so large that the convergence of 
Si (jcia) becomes too slow to justify the approximation made above ; and then th« 
formulae (5*4) to (5*6) provide the solution. It will be noticed that when K is 
complex, the equation 

cos e = A ( g^3 +2) M' 

will not usually give a real value for 6 : and so there is usually no direction in which 
the scattered wave is zero* (or of order (/ca) 7 ). 

It may be of interest to note here that experimental work on the scattering 
of light by fine particles has been carried out with silver particles suspended in 
water. f The corresponding values of tea seem to vary from ^ to 2, and the value of 
v/K is taken as 0*2 — i (3*6); thus the approximations in (ii.) are not sufficient to 
calculate either S w (kc&) or & n (K 1 a) with any accuracy. J In actual fact it proved 
necessary to use Lord Rayleigh's exact formulae, equivalent to (47) above, and to 
go as far as n = 4 in the series. § 

§ 6. Case of Large Perfectly Conducting Spheres. 

Before proceeding to the final formulae, it will be convenient to state certain results 
given by Macdonald|| for the values of the functions S n (z), E n (z), when both n and 
z are large. 

* For the case in which K - 1 is small this conclusion is given by G-. W. Walker, ' Quarterly Journal 
of Mathematics/ vol. 30, 1899, p. 217. The fornfulse given on that page agree with (5*6), when K- 1 is 
small; but the more general formulae on the preceding page do not agree with (5*6) completely. I have 
not succeeded in tracing the discrepancy on account of the fact that G. W. Walker has omitted some of 
the details of his preliminary calculations. 

t E. T. Paris, 'Phil. Mag./ vol. 30 (Ser. 6), 1915, p. 459. 

J To obtain an accuracy of 1 per cent, in Si (z) by retaining two terms of the series only, it must be 
supposed, that \z\ does not exceed 1 *3. 

§ E. T. Paris, be. ciL, p. 472. 

|| 'Phil. Trans.,' vol. 210, A, 1910, p. 134: the formulae are due to L. Lorenz originally. Avery 
interesting method of deriving the results is given by Debye ('Math. Annalen,' Bd. 67, 1909, p. 535). 

2 E 2 
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Provided that z—(n + jf) is of an order higher than z lj % the formulae are 



(61) 



< 



f E n (z) = -Re-'*, S„ (s) = B sin ^, 

where 

R 3 
and 

cos a 



l/sin a, 



\Js = 0sin a + |-7r— (n + Jr) a. 



= (^+i)M 



We shall need also the corresponding formula for S' B (z) and E' B (z) ; it will be seen 
that 



E' w (2) __ 1 e?R c?^ __ 1 ( x cos a da 



— L 



~E n (z) R dz dz WV sin 2 adz 



+ < > 



because (3*86) gives 



Hence 



/• 



(6-2) < 



where # 



R 2 ^ = 1. 



E n (2) 1 / , , \ 

-^ = __(tan x + .) 



£ 6 



_.i X 



R 2 cos x 



, _ 1 cos a c?a __ x cos 2 a 
tan x ~ 2 _-„ 2 ~ — 2 r ~ ir ~ 



V. 



sin 2 a <fe 



z sin a 



Similarly, we find that 



(6-3) 



S'. (z) _ 

8 n (z) R dz 



1 dR , , , d^ I/4 , . \ 1 cos(\^+x) 

BT-jT +cot ^ ^7 = ™ (-tanx + cot ^) = ^ • ., * 



dz R 



R 2 sin \fr cos x 



The formulae to be used finally are those for A n and C n , given in (471), thus 
we take 



■A-n — 



E'„ (s) 



G n 



E„ (z) ' 



where z now denotes kol. It follows from (6'1) above that 

C„ = -sin ire* 1 * = +ii (e 2 ""-l), 
and using (6'2) and (6'3) we see that 

» 

A n = -^ + * } cos(^+x) = -^{e 2t( ^ +x) + l}. 

It is now an easy matter to write down an approximation to the functions M and 
N defined in (4*8), provided that 6 is not near to or ?r. Under these conditions we 
can take the approximate value 



P n (cos 0) = /y/i 



2 



yiir sm Oj 



COs{(n + 4-)0-|-7r} 



* Under our conditions oc is not near to zero and z is large, so that x is small. 
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giving a corresponding approximation 

F (cos 6) = 2n+1 sin {( n +2-)fl-i^} 
" V ; sin0 v/(2n7r sin 0) 

Ihus 

sin , J- + 1 A ,p-, (cos „ . _ 1 , {e M.«. + i } (<! -_e-), 

where 

. = <»+*)*-*,, 

and we have replaced (n+^) 2 /{n (n+l)} by unity, because n is large. 
Thus, on putting x = 0, we find that (4*8) gives the approximation 

(6*4) M = -cos - — 2 —rnr — =~m (e** + 1) (^-e~ lw ). 

jct v / (2^7rsm0) 

Similarly, we get the formula 

(6*5) N = +sin0- — Z //o 6 . M (e at *-l)(e w -e" w ), 

v «r */(2mr sin 0) 

With series of this type, the leading part is found by making the index of the 
exponential stationary (regarded as a function of n). Now in both M and N there 
is one index only, a- = 2\^ + nw— w, which can be stationary : and the condition is 

an an 
Now, from (6*1) 

d\lr r / t 1U da 

-^ = {zcos a -(n +3r )}^-« = -a, 

and so the leading terms in (6*4) and (6*5) arise from taking 

2a = 7T— 0, or n + % = z sin |-0. 
The corresponding value of the index or is then 

cr = 2^+n7r— to = 2z sin a+|-7r— (2n + l)a + mr— (n + ^fl+iw, 

To determine the form of the index near to this special value of n we take 

a <j q eta Z Jj 



dn 2 dn z sin a z cos ^6 

and then we find the approximate formulae 

f , (n-n ) 2 ' 

<r = Vol fT> 

(6*51) < , 2fC0S ^ 

1 7 where 

^o + i" = 2sin|-0. 
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The leading parts of M, N are accordingly given by the approximations 



(6*6) 
where 



e 



•IKf 



e 



— mr 



M = + cos </> — - Q, N = —sin <p Q, 



kT 



kT 



Q = 2 



e 



tar 



x/(2u7r sin 9) 
and a- has the value given in (6 '51). 

The value of Q is approximately equal to the integral 



»oo 



(6'61) 
where 



oo 



e l<r dx 
*/(2n Q Tr sin 9) 



n = n +x y or = or ~\-x 2 j(z cos %0). 



Thus, approximately 



Q 



e 



L<T 



\/(2n 7r sin 6) 



VlirZ cos 



#' 



e 



i(<r Q + irr) 



ie 



t (2z cos f 0) 



2 sin 4$ 2 sin A-0 



2 



2 



Accordingly, to the same degree of approximation, we can take 



(6-62) 



de 



X. l yo^ l Z COS id 



Then the components of force are given, as in (4*9), by 



r 



(6*63) 



< 



X ■ = + Cy == 



Z = -c/3 = 



1 SN 



3M 



30 sin 6 d(p 

. 1 3M3N 

+ 



e 



— t«r 



COS — — -r~ , 



, — IKf 



V 



sin 9 d<p d9 



, e~~' dQ 
sin <h ~ 9 

7 kt d9 



where differential coefficients with respect to <p are small compared with those with 
respect to 0, and so have been rejected. Thus, using (6*62), we have the approxi- 
mations to the forces in the scattered waves 



r 



(67) 



Y = + cy — cos <p — er > 



2r 



-< 



Z = -c/8 



V. 



sin 6 — e 2lKacmi6 ~ lKr 
Y 2r 



assuming that 9 is neither near to nor to it. 

When 9 is small, the approximation to P n (cos 9) must be taken as 



and so 



P n (cos 9) == J {(2n + 1) sin -J0}, 
sin £P n (cos 0) = (n + i) cos |-0 J x {(2w+ 1) sin ^0}, 
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Thus we now write 



sin 6 



2n+l 
n (n+1) 



A n P' n (cos 6) = -i (e 2t * + 1) cos |0 J x {(2n + l).sin £0} 



and, proceeding as before, we are led to the conclusion that a is near to \ir ; thus the 
value of n/z is small, in the parts of the series which contribute the principal part of 
the sum. Then we can replace (6*1) by the approximate formulae 



;. 



(6*8) < and 



Thus 






Hence the approximation corresponding to (6*4) is now 



(6*81) 



M = -* cos £0 cos <j> - — S e^+'fc+^'J! {2n+ 1) sin £0}. 



AT 



In like manner the value of N is found to differ from (6*84.) only in having + sin <f> 
as a factor instead of — cos <p. 

In the series (6*8 1) the value of n may be supposed to vary from to oo ; and so 
we obtain the principal part of the sum by using the integral 



(6*82) 



/•CD 

e 21 * e^J, (2£ sin |0) d£ ; 

Jo 



and when sin J0 is very small the value of (6*82) is approximately equal to 



Thus 



e 2lz (%tz sin |-0). 



/* 



(6*83) 



M = +^z sin cos <pe 



2tz 



e 



•IK?* 



KT 



< 



and 



N = — \z sin sin <j>e 2lz 



e 



• «r — 



v 



KT 



Accordingly the components of force are now found to be 



r 



(6-9) 



< 



Y= + Cy 



Z] = -c/3 



3M 1 3N 
30 sin 6 3^ 



1 3M , 3N 

+ 



1 S 

f Z COS <£ - 



2iz—f-<r 



KT 



V. 



sin 3^> 30 



i • J e 
■%z sin <£ - 



2t2-— t«r 



icT 
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These results in (6*9) are precisely the same as would be found by writing 6 = 
in the approximations (6*7), and accordingly the formulce (67) do remain valid right 
up to the axis 6 = 0. 

When 6 is nearly equal to x, the calculation on the present lines becomes more 
difficulty and we shall accordingly obtain the corresponding approximation by a 
different process in the next section (§ 7). 

It appears that for the special value tea = 10, the formulae (6 '7) do give the forces 
with a fair degree of accuracy up to an angle 6 = § it ; the approximation in fact 
appears to be better than might have been expected. [See p. 176 above.] 

§ 7. Alternative Method, applicable to any Conductor whose Dimensions 

are Large compared with the Wave-length. 

It follows at once from Green's theorem that if u> v are solutions of the equations 

A 2 u+k 2 u = 0, A 2 v + k 2 v = 0, 

at points within a closed simple surface S, then 

(7-d ;(„£_,£)*_* 

where the integral is taken over the surface S (dv being the element of outward 
normal), and it is supposed that u 9 v are both free from singularities in the interior 
of S. 

Similarly if u has no singularities and v is a solution which behaves like e~ l * R /Il 
near a particular point P (R denoting the distance measured from P), we see that 

(7*11) \\ u ^ v —) dS = — 4:7ru P , 

J \ cv dv I 

provided that P is inside the surface S. 

Equations of similar forms apply when the space considered is outside the surface 
S ; but then the sign of the last equation (7*11) is reversed, giving 

(7-12) M u to- v ^\dS = +4tt^ p ; 

it is then necessary to assume also that at infinity, u> v both correspond to 
divergent waves (unless it is known that u tends to zero more rapidly than l/r). 

* Compare Macdonald, loc. cit., pp. 120-122. The cause of the difficulty is to be found in the fact 
that now the stationary value may be expected to arise from values of n for which a is small. Then n is 
nearly equal to #, and in all such cases more complicated analysis is inevitable. In fact the approximations 
to S n (#) and E n (z) require to be modified by different formulae corresponding to the cases n > z, n < z 
and to the cases in which j n ~ z | is of order z 1 /* or of lower order. 
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Let us now consider the problem of waves incident from some source (or sources) 
and reflected from the surface S. Let u denote any Cartesian component of force in 
the incident wave, and let the point P be outside S. Then (7*1) applies, and so 



i 



u 



dv du 



dv 



v 7T- )cB3 

ov J 



0, 



because u has no singularity inside S. 

If v! denotes the corresponding component of force in the reflected wave, we have 
from (712) 



dv/ 



because u f has no singularity outside S (and u' will correspond to a divergent wave 
at infinity). 

By addition we have the result 



(713) 



4<7rt/ P = 



ov dv ' '' 



where v is taken to be e~ tKR /E,. 

Now u + vf = w gives the corresponding component of force in the complete wave ; 
and this accordingly satisfies certain known relations at the surface S (the exact 
form depending on the physical properties of S). It must, however, be clearly 

understood that we cannot usually obtain both w and -5— by any simple methods, any 

more than the analogous problems of electrostatics can be solved by a mere appeal 
to Green's Theorem. 

However, we can obtain an approximate solution, suitable to the problem of short 
wave-lengths, by assuming that near the reflecting surface, the character of u f can be 
determined from that of u by the rules of elementary geometrical optics. Thus we 
treat the reflected wave as derived from the incident by simple reflexion in the 
tangent-plane at the point of incidence. Making this hypothesis it is an easy matter 

to construct both w and -~ when the form of u is given. 

It will be noticed that we shall have w = 0, -77- = at all points within the 

ov 

geometrical shadow ; and so the final integral (7*13) extends only over the illuminated 

side of the surface S. 

Suppose now that we consider electric waves incident on a simple convex 

conducting surface ; and take an origin on the surface such that OP is the 

reflected ray (in the sense of geometrical optics). Take the plane of incidence as the 

plane of yz, and the normal at as the axis of z. 

YOL, CCXX. — A. 2 F 
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Then in the immediate neighbourhood of O we can represent the incident wave by 
the components of electric force 



(A, B, C)e 



ik (my + nz) 




where m 2 + n 2 — 1, so that m is the cosine and n the sine of the angle of incidence. 
We have also the relation 

Bm + Oi = 



because the electric force is perpendicular to the incident ray. 
The corresponding reflected wave has the components 



lk (my—nz) 



where 



and 



so that 



(A', B', C) e 



A + A' = 0, B + B' = 0, 
B'm-C'n = 0, 



C ; -C = 0; 



these results follow by making the tangential components of force zero on the 
tangent -plane z = (instead of at the surface). 

It is now clear that, at the point O, the components of the total force will be 



(7'2) 



(0, 0, 20), 



and the normal differential coefficients of the total force will be 



(7*21) 



lku(A-A\ B-B ; , C-C / ) = 2actc(A, B, 0). 



dw • 



We shall now insert these values for w and — in the general formula (7*13), 

including also the factor e LK{mj/+nz \ on account of phase-differences at points near to O, 
Since reflexion takes place only from the immediate neighbourhood of O, the error 
introduced by this simplification will be small 
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Since the co-ordinates of P are (0, —mr, nr) y the value of B, is given by 

B 3 = x 2 + (y + mr) 2 -\- (z—nr) 2 , 
= r 2 + 2r(my— nz) + x 2 +y 2 + z 2 , 

where (x, y, z) is a point on the surface near to O. Thus, when r is very large in 
comparison with the dimensions of the surface (as we assumed in the previous 
investigations, §§ 4-6), we can use the approximate formula 

(7*3) R = r+my—nz. 

Thus we can write in (7*13) 

p — iK(r+my—nz) 

v = = V{i e- lK{my - nz \ 

r 

^\ 

if v is the value of v at O. Then the most important term in — is seen to be 

cv 



dv 
Zz 



— = ucnv ( /r ut{m *-n*\ 



Accordingly the components of electric force in the reflected wave will be given by 
the approximation 



r 



(7-4) 



X = -4* 



:7T 



(2cKfiA) e 2tKnz dS, 



< 



4x J 
Z = + -^ f(2«nC) e 2w "<2S. 

47T J 



To evaluate the integrals in (7*4) we must write out the equation to the surface in 
the approximate form 

2z = —(ax 2 + 2/3xy + yy 2 ). 
Then* 

J lkU <sj(ay— /3 2 ) 

Now ay--/3 2 is the absolute (or Gaussian) curvature of the surface at the point O ; 
and since the surface is supposed convex, we represent this curvature by \\p. 

* This is most easily found by taking the integral as lim L-^+^w^+tf^+Y^dS. 

2 F 2 
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Thus we get 
(7'41) iKn\e 2tKnz dS = irp. 

Hence, using (7*4), we see that the principal part of the reflected wave is 
given by 

(7-5) (X, Y, Z) = (-A, -B, +C) f r e-". 

In order to interpret (7*5) for any axes of co-ordinates, we need only notice that 
(—A, — B, -fC) represents a force numerically equal to the force in the incident 
wave; and that the new force is perpendicular to the reflected ray, arranged in such 
a way that the tangential components are opposite to those in the incident wave. 

We can apply the formula (7*5) to the problem of § 6 at once ; clearly p = a. 

The point of incidence corresponding to the scattered wave (#, <}>) is given by (|r0, ^). 



X 




Then the incident wave at O has the components of electric force 

— cos 0e tKa cos ^ e in the plane of incidence, 

and 

+ sin <f>e lKa cos ¥ perpendicular to the plane of incidence. 

Further, the r of formula (7*5) is measured from O ; to compare with § 6, we take 
r to be the distance CP, measured from the centre of the sphere. Thus we are to 
replace r in (7*5) by i — acosJk?. Accordingly the components of force at P, in the 
reflected wave, are 



a 



and 



H cos 0e 3MtaoOB * fl " uer perpendicular to r in the plane ZCP. 

2r 



nln A. ^>2l(ca COS §0 — IKT 



sin ^e 21 ^ 008 ^- 1 ^ perpendicular to the plane ZCP. 



2r 
These results agree with (67) and (6*9) of § 6 above. 
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It is easy to modify the general formulae (7 ' 5) so as to cover the case of waves incident from a point 
source (say at distance s from the point of incidence). 



X 




Then 



Sl 2 = x 2 + (ms - y)*<+ (ns - zf 
= s 9 - 2s (my + nz) + x 2 + y 2 + z 2 



Thus with the usual approximation of geometrical optics 



Similarly 



si = s-(my + nz)+ — {x 2 + (ny-mzf}. 

AS 



R = r + m y-nz+ — {a? + (ny + mz) 2 } . 



Now z is of the second order in comparison with x, y ; and so we can write 



Thus here we find 



where 



Si + R = $ + r - n(ax* + 2fixy + yy 1 ) + ( —• + — ) (x 2 + nY). 



J iKTur 



cr 



1(1 -IV- 1 



w 



r h r 2 being the distances of the focal lines (of geometrical optics) from the point of incidence. 
Thus now the principal parts of the reflected wave are given by 



(-A, -B, +C)e~ 
assuming that r is not close either to ri, or to r 2 . 



-/{(' -S('-0" ! 



The results of the foregoing analysis depend on the tacit assumption that n is not 
zero ; and as a consequence the character of the approximations will change when n 
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is small. That is, near the edge of the shadow, in the ordinary phrase of geometrical 
optics. 

In the application to the sphere, the region excluded by this condition corresponds 
to values of 6 nearly equal to 7r ; and in this region the specification of the scattered 
wave by means of Cartesian co-ordinates seems simplest. 

We consider then the incident wave as specified in Cartesian form by 

(7*6) X = -e l % Y = 0, Z = 0, 

and consider the approximation to the reflected wave incident at the point (al, am, an) 
on the sphere for which the direction-cosines of the normal are I, m, n. The expressions 
will be of the form 

(7*61) (X', Y', 71) = (A', B', C)e wC 

where (treating the tangent-plane as the reflecting surface) 

f = z+p (a—lx—my—nz) 
and p is determined by 

( — lp) 2 +(—mp) 2 + (l—np) 2 = 1. 
Thus 

p = 2n. 

Further, the resultant of (X, Y, Z) and (X' ? Y ; , 71) at the point of incidence must 

be along the normal ; and so 

A 7 -! = B' = <7. 

I m n 

Also (A', B', C) is perpendicular to the reflected ray ; and so 

A' (-2nl) + B f (-2nm) + C (l-2n 2 ) = 0. 
Hence 

A 7 -! B' C 2nl = ___ 2/ 



(7*62) < 



I m n n—2n 

or 

A' = 1 - 2Z 2 , B' = - 2lm, C = - 2ln. 



The components X + X', Y + Y ; , Z-fZ' at the point of incidence are accordingly 
equal to 

(7*63) -2(l\lm, ln)e LKan . 

We have still to evaluate the normal differential coefficients, which are found to be 
(7*64) ^(-1-A', -B', -C')e lKan = 2um(f , -l, Im, ln)e LKan . 

The value of B is now seen to be given by 

B 2 = (x— al) 2 + (y— am) 2j r(z— an) 2 = r 2 — 2a (Ix-h my -\-nz) -{- a 2 . 
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Thus, when we regard a/r as small, we may take 



(7'65) < 



r R = r— a (I sin. #cos <p + m sin 6 sin <p+n cos #), 

and 



3R 

— = — (I sin cos (p+m sin sin c/>+n cos 0). 

01/ 



3R 

Now, in applying (7*13), — occurs only in the coefficients and not in the exponential 

ov 

3R 

index : thus we can get the first approximation by putting 6 = -k in — ; this gives 

ov 

the value 

(7*66) ~ = n. 

ov 

Thus, to our degree of accuracy 

47T \ di/ ov/ 47r\ ov di/ 

V I , dw" 

IKUW+ — 
OV 



:7T 



We can now substitute for w and — the values given by (7*63) and (7*64) : it will 

ov 

be seen that the components parallel to y, % give zero (to this order), and that the 

component parallel to x gives 

iku e iKian ~ B) 



2-7T r 



Accordingly the reflected wave is given by 
(7*67) X = — [ndSe tKian ~ n \ 

where R is found from (7*65) and the integral extends over the positive hemisphere. 
For the purpose of integration we write 

I = sin ff cos <£', m = sin 6 f sin <f>\ n = cos 6'. 

Then in (7 '67) we have 



an- 



.j{ — —r + a {(l + cos #)cos fl'-fsin #sin O f cos (<f> — <p')}. 



The integration of (7*67) with respect to </>' can be carried out at once, because 

/*27T 

(7'68) e «orin^Binroo8(*-^) ^/ = 27r J ( Ka s in s i„ #'). 

Jo 
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Our result accordingly becomes 



(7'69) 



X 



LKOL 



T 



■ tar 



2 Va(i + cos0)cos0'j o (ica sin 6 sin 6') sin tf cos 6' dtf. 





An exact evaluation of the integral would be troublesome ; but it can be trans- 
formed by integration by parts. This process leads to a series of which the first 
term is 



(77) 



A. — e . 

r ica sm 6 



The formula (7 '7) will represent (7*69) sufficiently accurately if (l + cos #)/sin 9 = cot ^6 
is regarded as small ; and this is the correct assumption here, since 9 is supposed to 
be nearly equal to ir. ' 

In the special cases ko, = 9, 10 it appears that the approximation (7 '7) represents 
the scattered wave sufficiently well within a cone extending to about 10° from the 



axis* 



